Nonlinear ion-acoustic (lA) waves driven in a 
cylindrically symmetric flow 
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Abstract By employing a self-similar, two-fluid MHD 
model in a cylindrical geometry, we study the features 
of nonlinear ion-acoustic (lA) waves which propagate 
in the direction of external magnetic field lines in space 
plasmas. Numerical calculations not only expose the 
well-known three shapes of nonlinear structures (sinu- 
soidal, sawtooth, and spiky or bipolar) which arc ob- 
served by numerous satellites and simulated by mod- 
els in a Cartesian geometry, but also illustrate new 
results, such as, two reversely propagating nonlinear 
waves, density dips and humps, diverging and converg- 
ing electric shocks, etc. A case study on Cluster satellite 
data is also introduced. 

Keywords nonlinear waves; ion-acoustic (lA); MHD; 
satellite 



1 Introduction 



Nonlinear plasma theory and approa ches have been de- 



Then. nonlinear surface waves were found to shift the 



veloped for more than half a century (jSagdeev fc Galeev 



196a iDavidsonI 119721) . Important problems, such as 



the excitation, propagation, and effects of nonlin- 
ear waves, have been extensiv e ly stu died since 1970s, 
see, e.g., Ilnfeld fc Rowlandsl (|2000l ). One branch 



lies in cylindrically symmetric plasma systems which 
are ubiquitously observed in geo-space by a multi- 



aeers ( Pimenta et al. 


2001). rockets (Earle et al. 


1989 


Moore et al. 1996). and satellites ( 


Pickett et al. 


2004 


Vaivads et al. 


2004; 


De Kevser et al. 


2005). 



In the cylindrical frame where B is along axial z- 
direction, small- amplitude ion acoustic f lA) waves wer e 
firstly studied (jMaxon fc Viecelil Il974 iMaxonI Il976l) . 
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radius of a plasma cylinder (iGradov fc Stenfid 11983 



Gradov et al.l 119841 ). After that, the deformed ampli- 
tudes and propagation velocities were confirmed to be 
different from those i n the semi-infinite limit of plasmas 
( Gradov et al.lll985L Il986[ ). In addition, the equations 



governing evolution of surface and body waves wer e 
obtained |Molotovshchikov fc M. S. RudermanI [l987l) . 
Furthermore, in unmagnetized systems, a set of non- 
linear equations was solved to describe the temporal 
change of t he electron density in strongly no nlinear sur- 
face waves (|Stenflolll99dlYu fc Stenflolll99lh . The work 
was followed by a generalized study including the ro- 
tation effects of a t i me-de pendent rigid plasma body 
( Stenflo fc Yul Il992l . Il995l) . By contrast, in magne- 



tized systems, nonlinear lA parallel-propag ating elec 



trie field st ructures were recently calculated (jShi et al 



20011 120051) , and bipolar electric field structures were 
verified to be able to originate froni either lA or ion 
cyclotron (IC) waves |Shi et al.ll2008h . 

By examining observations by high-resolution satel- 
lites (e.g.. Wind, FAST, Polar, Cluster), we found 
that, in addition to the bipolar shapes of nonlinear 
structures, there exist two other well-known nonlin- 
ear electric field envelops: sinusoidal and sawtooth. 
These three shapes were first recorded by S3-3 in 

Since then, they 
satellites such as S3-3 



the 1970s (jTemerin et al.l 119791 ). 
been detected 



numerous 



(Temerin et al. 19821). Viking i 


Bostrom et al. 19881). 


Geotail (Matsumoto et al. 1994 


). Wind (Bale et al. 


1998 


). FAST (Ergun et al. 19981 McFadden et al. 


1999, 


2003 


). Polar (Mozer et al. 1997; 


Franz et al. 


1998; 


Bounds et al. 


19991 Cattell et al. 


1999; Franz et al. 


2000f). and Cluster ( 


Pickett et al. 2004. 2005). 



The formation of these three nonlinear structures 
have been stu died extensively in Ca rtesian coordinates. 
For example, iTemerin et al.l (|l979[ ) firstly reproduced 
these shapes by solving a set of fiuid equations in the 



2 



ion-cyclotron (IC) / ion-acousti c (lA) regime. In a 



unified work, Lee fc Kan ( 198l[) obtained these two 



important nonlinear electrostatic waves (it is worth 
to mention here that the authors also mentioned a 
third type of the nonlinear waves: "ion-acoustic soli- 
tons" . By redoing the calculations, we can easily 
see that it is not a new type but the simple waves 
with a longer period for p eaks to occur). T heir stud- 



ies were followed by, e.g., iNakamura fc Sugai 



Chatteriee fc Rovchoudhurvl (119971): 



Jovanovic fc Shukla (200C) 



Partic u larly. iReddv et al 



Das et al 



Mamun fc Shukla 



1996 ): 



2000 ): 



20021) 



(|2002l) : iBharuthram et al. 



I 2OO2I ): iMa fc Hirosel |2009l) not only verified that the 



nonlinear structures originate from a coupling between 
the IC and lA modes, but also obtained the thre e wave- 
forms reported firstly by iTemerin et al.l (|l979l ). The 
authors found that the nonlinear shapes change from 
a sinusoidal IC/IA structure of small amplitudes at 
low Mach numbers to sawtooth and then bipolar/spiky 
waveforms of large amplitudes at higher Mach numbers 
in parallel propagations. 

In cylindrical coordinates, however, there was, to our 
knowledge, only one study on the contribution of cylin- 
drically symmetric plasmas to the emergence and prop- 
agation of nonlinear str uctures measured by sate llites. 
This study was done bv lTrines et~al] |2Q05l [20061) . The 
authors offered a matched picture between the Cluster 
data of nonlinear waves and the kinetic modeling of the 
wave excitation in the perpendicular plane of the local 
magnetic field B. Inspired by their work, we pay at- 
tention to parallel-propagating nonlinear waves, aiming 
at clarifying if or not the well-known three envelops of 
nonlinear structures can be driven in cylindrically sym- 
metric flows, and if possible, obtaining new features of 
nonlinear waves in such a system, such as density holes 
and humps. The organization of the article is as follows. 
Section 2 set up a two-fluid model by employing a set 
of self-similar MHD equations. Section 3 introduces a 
parameterized analysis of nonlinear waves by numeri- 
cal calculations. Several features of them are obtained. 
Section 4 compares the results with that in a Carte- 
sian frame, and have a case study on the Cluster data. 
Section 5 presents a conclusion. 



2 Two-fluid model 

In order to provide the most basic picture for the prop- 
agation of electrostatic nonlinear waves driven in cylin- 
drically symmetric plasma regions, and be able to il- 
lustrate clearly the modulation of nonlinear structures 
by various input parameters and boundary conditions, 
we focus on FAST/Cluster orbits where the plasma /3 



is much less than 1 (such as, the auroral acceleration 
regions, bow shock, magnetopause), and employ a two- 
fluid model in a cylindrical frame (r, z) with B^Se^ 
(where B is constant, is the unit vector along z axis). 

The model takes into account isothermal electron 
and ion fluids, with vtc^vti, where WTe= \/2fcs7e7"T^ 
and VTi=\/'2kBTi/mi are the electron and ion thermal 
speeds, respectively, in which fc^ is Boltzmann's con- 
stant, Tg and Ti are their temperatures (assumed con- 
stant), respectively, and TOe and are their masses, 
respectively. In this study, we neglect the electron iner- 
tia because of me<C rrii. Both the momentum equation 
and the isothermal equation of state of the electron fluid 
provide Ne=e^ , where N^~ni./nQ is the dimensionless 
electron density in which rie and rig are the dimensional 
electron density and the uniform background plasma 
density, respectively, and <^=eLp / (kBT^) where ip is the 
perturbed electrostatic potential. 

The ion fluid is described by a set of MHD equations 
expressing ion continuity, momentum (containing the 
equation of state), and the quasi- neutrality: 
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in which Ni=ni/nQ (where n.i is the ion density), 
T=Qit (where fli=eB /rrii is the ion gyro-frequency, 
and t is time), Vr=Ur/cs,V^=Uif,/cs,Vz=Uz/cs (where 
{ur,U0,Uz} is the ion velocity components, and Cs ~ 
yJkBTe/mi is the ion acoustic speed), R=r/pi, Z=z/ pi 
(where pi=Cs/^i is the ion gyro-radius), and C = 
1 -|- i;|,j/(2Cs). Note that (1) due to the symmetric na- 
ture of the cylindrical system, all derivative terms along 
the 0-direction do not occur; and (2) the parameter R 
represents the radius of curvature of the flow stream- 
line intersecting the magnetic field lines on which the 
equation is going to be solved. 

In this system, linear lA and IC waves can be excited. 
By linearizing Eq.(IT]), we obtain a;^=^fc^c^ (where k is 
the amplitude of the wave vector k) in the parallel di- 
rection and a;|=f2f in the transverse plane. Superim- 
posing upon these background oscillations, there exist 
nonlinear waves the features of which can be obtained 
by solving a set of self-similar equatio ns of Eg. (HI) via in- 



troducing a sel f-similar parameter X (|Lee fc Karj|l981 
Shi et ahlkofil : 



X^^R+^Z- 
M M 



(2) 



where M is the Mach number, Q;i=sin6', and q;2=cos6' 
in which 6 is the inclination angle between the propaga- 
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Fig. 1 (Color online) Ion density Ni and wave-field strength Ei of sinusoidal nonlinear waves excited under V|/>o=0.05 at 
M=2 and i?=l. Upper panel: Fast nonlinear wave; Lower panel: Slow nonlinear wave. Inserted in each panel is the speed 
of the simple wave propagating along B. 



tion direction and the magnetic field. Using self-similar 
transformations, i.e., 



d_ 



d 



d_ 
dR 



ai d 
MdX' 



dZ 



a2 d 
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(3) 



We concentrates on parallel-propagating nonlinear 
waves, i.e., 9=0. This means ai=0, and Q!2=1- There- 
fore, we have 



_d_ _a_ _ 1 d 

dR~ ' dZ ~ MdX 



(4) 



the first expression of which indicates that R is inde- 
pendent of X. We thus obtain a set of four self-similar 
equations of nonlinear waves propagating along B: 
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(5) 



in which R behaves only as an input parameter to 
represent the existence of geometrical effects (namely, 
the centrifugal and Coriolis forces) in curvilinear flows. 
Note that these effects are absent in rectilinear systems. 



Using boundary conditions $|a'=o = 0, the third 
equation of Eq. ([5]) can be integrated to give 



(6) 



where 5'=2C$/Af^ and Vz|x=o have two initial values: 
and 2 A/. This leads to a new set of equations as 
follows: 



^ dX \R 



dV^ 
dX 
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M dX 
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where Ei is the dimensionless electric field amplitude 
of the simple wave with a unit of Eq=CsB. 



3 Parameterized analysis 

3.1 Prerequisite of the excitation of nonlinear waves 

Cylindrically-symmetric fiows are azimuthal with V^o 7^ 
at X^O. Other parameters satisfy equilibrium con- 
ditions, i.e., V^o='i'o=0. A nonzero V^q is important. 
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Fig. 2 (Color online) The same as Fig[l]but sawtooth nonlinear waves excited under V'</,o=0.5. 



This can be seen from the first two equations of Eq.Q. 
They provide 



VrdVr + V^dVA, = 



or, 



yI\x=o 



(8) 



(9) 



Clearly, if V(/,o=0, both Vr and are zeros at any X, 
and thus, the last equation of Eq.([7]) gives ^"=0 at any 
X. In this case, no nonlinear waves can develop. 

As a result, only under nonzero V^o conditions is it 
possible to trigger nonlinear processes. Luckily, this 
condition is always met naturally in toroidal flows in 
geospace. 

3.2 Evolution of nonlinear waves 

At three different levels of V0o=O.O5 (weak), 0.5 
(medium), 0.95 (strong), respectively, FigsllJS] illus- 
trate, respectively, the three well-known structures of 
the nonlinear wave-field Ei^ namely, sinusoidal, saw- 
tooth and spiky or bipolar, as well as the accompanied 
nonlinear- wave density iV^, under il/=2 and R=\. Ev- 
ery figure contains two panels, each of which is inserted 
into a Vz-curve to express the propagating speed of a 



nonlinear wave along the axial direction. Note that 
the fast wave is always propagating in the direction 
of B (Vz>0), while the other in the opposite direction 
(K<0). 

The three figures illustrate that the weak fiow drives 
sinusoidal nonlinear structures with small amplitudes 
(l-Eilmaa; =0.03), the medium flow drives sawtooth 
structures with medium amplitudes (|-E'i|maa:=0.3), 
and the strong flow drives spiky (bipolar) nonlinear 
waves with high amplitudes (|£'i|maa:=0.5)^ This re- 



sult r e produces the results obta ined by iReddv etal 



( 20021) : iBharuthram et al.l m Cartesian coordi- 

nates. This indicates that the features of nonlinear 
waves are strengthened increasingly with stronger drifts 
of azimuthal flows. Especially, in the last case, the two 
panels in FigO reveal diverging shocks^® (a negative 
electric fleld followed by a positive one) of the fast 
nonlinear-wave packet in the upper panel, and con- 
verging shocks^"^ (a positive electric fleld followed by 
a negative one) of the slow packet in the lower panel. 
In addition, more calculations with a changing M ex- 
pose that under M>1. the nonlinear- wave density Ni 
is never larger than 1, but goes to a minimum which 
is smaller at a faster drift, meaning density holes are 
formed and their boundaries result in nonlinear waves. 
For example, the hole is 0.876 in amplitude with a 



5 




smaller \Ei\max at V^o=0.05, while it becomes 0.015 
with a larger l-Bilmax at V^o=0.95. We note that V^o 
cannot be larger than 1, that is, the azimuthal flow 
speed is unable to surpass the local acoustic speed. 

3.3 Influence of R. 

In a cylindrical system, the frame effects produced by 
the centrifugal and Coriolis terms in the momentum 
equations decrease with radius. They should have a 
direct influence on the structure of nonlinear waves ex- 
cited in the azimuthal flows. Fig|4] expose the role 
played by R. At a larger radius R=5 compared to 
FigEl the nonlinear waves change their appearances 
from a spiky (bipolar) shape to a sawtooth one, with a 
lower amplitude \Ei\max=0-^'2 but a higher frequency. 
At another radius, i?=10, a calculation shows that si- 
nusoidal shape emerges, with a very small amplitude 
of \Ei\max~0-OQ- This indicates that the features of 
nonlinear waves arc weakened when going farther from 
the center of the cylindrically symmetric flow. It is 
thus reasonable to propose that it is easier to detect 
nonlinear waves closer to the symmetric axis where the 
influence of the centrifugal and Coriolis forces are more 
dominant. 



3.4 Criterion for density holes and humps. 

The Mach number cannot always be larger than 1,"^^ 
referring to that M<1 is also possible. Fig (S] illustrates 
a case under conditions of FigEl but at A/=0.5. In- 
terestingly, different from the density holes, there are 
now density humps coming into being: the nonlinear- 
wave density iV^ is never smaller than 1, and goes to 
a maxmum 1.133. In addition, the bipolar nonlinear 
wave becomes denser with a much higher frequency of 
3, than that of 0.04 in FiglS] However, the amplitude 
of Ei decreases from 0.5 to 0.35, meaning the nonlinear 
feature is weakened. Other calculations confirm that 
density holes occur at M>1, while density humps ap- 
pear at M<1. 

4 A case study 

We have chosen a cylindrical geometry to study the 
excitation of nonlinear structures modulated by differ- 
ent values of the initial azimuthal flow speeds, the dis- 
tance to the symmetric axis, and the Mach number. 
We not only reproduced the three familiar shapes (i.e., 
sinusoidal, sawtooth, and spiky or bipolar) of nonlinear 
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Fig. 4 (Color online) The same as Fig[3]but sawtooth nonlinear waves excited at a distant radial position _R=5. 



waves observed by satellites, but also found that nonlin- 
ear waves supported by cylindrically symmetric plasma 
flows have both fast and slow branches, both converging 
and diverging electric shocks, and, both density humps 
and dips. This is different from the rectangular flows. 
By solving Eq. ^ in a Cartesian frame, we obtain that 
the w ave potential $ satisfies [c.f., e.g., iMa fc Hirose 
( 20091) ] 



e — 



%/2 



r)+V[(^ 



-12- 



(10) 



where ^ is the ratio between ion and electron initial 
temperatures. Taking ^ = 0.1 and M=l. 14, 1.16, 1.28, 
respectively, produces a result shown in Figl6] Though 
the figure keeps exhibiting the three familiar nonlinear 
structures, only diverging shocks are seen to be driven. 

The cylindrical model and calculations presented in 
last Sections are useful in gaining physical insights into 
the excitation and propagation of observed nonlinear 
waves. Figl7]is an example, exposing snapshots of non- 
linear waves within the plasma sheet on auroral field 
lines, taken by the four Cluster satellites when they 
encountered the high altitude auroral zone at a radial 



distance of -4.2 Re at 21:30MLT and ~ 62° ILAT 
(sout h) on 31 March 2001 at -06:44 UT (jCattell et al.l 
2003h . The panels show both converging shocks (SCI 



panel) and diverging shocks (SC2-4 panels), sawtooth 
structures (SC2 panel), and more complicated shapes 
of background hiss, at the four locations separated by 
— 1000 km and an —0.03 s interval from each satellite. 
The detected amplitudes of waves can reach over 750 
mV/m, the largest values ever reported in the outer 
magnetosphere. 

At that time, the mission was experiencing a strong 
magnetic storm. The magnetosphere was highly com- 
pressed djae_to_a2jntense injection of electrons —14 min 
earlier (jBaker et al.1120021 ). which could enhance the lo- 
cal magnetic field from a few mC of the background 



strength to tens even hundreds of mG (|Alfaro et al. 



2004 iTiulin et al] booi iHall et al.l boOOl) . Although 



the Cluster has no single probe measurements to give 
nonlinear wave speeds, all four panels reveal that they 
excellently record ed all the bursts by double probes of 



high resolutions (jGustafsson et al.l Il988t iCattell et al 
2003f ). The pulses are —2.5 ms or its multiples in 



time, in good agreement with our simulation value of 
X=10-30 by using the local orbital data of the Clus- 
ter. More important, SCI recorded an oppositely prop- 
agating nonlinear wave of converging shocks, different 
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Fig. 5 (Color online) The same as Fig[3]but iiiglier- frequency spilcy (or bipolar) nonlinear waves excited by a low Macli- 
number Ad— 0.5. Nonlinear waves have density humps. 



from other satellites which displayed diverging shocks, 
but at the same time. This is in good agreement with 
the illustration of two oppositely propagating nonlinear 
waves existing in a cylindrical system simultaneously, as 
shown by FiglS) This fact notwithstanding, we are un- 
able to find the wave speed data to double-check Vz of 
the wave detected by SCI, which should have a smaller 
amplitude than those measured by other satellites. 



5 Conclusion and discussion 

By using a two-fluid self-similar MHD model, we stud- 
ied features of parallel-propagating nonlinear waves 
driven in a cylindrically symmetric system. We not only 
reproduced the three salient shapes (i.e., sinusoidal, 
sawtooth, and spiky or bipolar) of parallel-propagating, 
lA/IC nonlinear waves, but also found following new re- 
sults: (1) the prerequisite to trigger the nonlinear waves 
is the nonzero azimuthal speed; (2) there are always two 
types of nonlinear waves: a fast one propagating along 
B (diverging shocks) and a slow one against B (con- 
verging shocks); (3) the distance from the symmetric 
axis influences the nonlinear features: the closer to the 
axis, the more pronounced the characteristics, and vice 



versa; (4) there exist density holes for M>1 and humps 
for M<1, the boundaries of which constitute nonlinear 
structures. 

This study is the first step to investigate the mech- 
anism of nonlinear wave-particle interactions and their 
effects on some unusual observed phenomena, e.g., 
transverse ion heating, broadband noi se emission, and 



magn etic holes (or bubbles, decreases (jTsurutani et al 
200,^ . It offers an alternative to expl ain inverted-!^ 



struc tures in beam-precipitating regions ( Pottelette fc Berthomier^ 
20091) . as to be contributed in a companion paper. This 



is based on the considerations as follows. 

When a field-aligned current is enclosed in the iono- 
sphere by the Pederscn current, plasma turbulence en- 
hances abnormal resistance locally which may bring 
about an electric field perpendicular to the magnetic 
field lines. This field drives naturally an azimuthal fiow, 
and thus provides, on one hand, a boundary condition 
to Eq.([T|). On the other hand, there also exists an- 
other magnetospheric boundary condition. Therefore, 
unlike the treatment employed in the present paper, the 
solution that yields a non-vanishing parallel potential 
drop has to be sought in a system where Eq.([T]) should 
be solved two-point boundary condition problem, 
rather than an initial condition one as was the case 
here. From the result of this paper, we know that the 
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Fig. 6 Parallel-propagating nonlinear waves developed 
from a rectangular flow. Although there still exist three 
types of shapes (sinusoidal, sawtooth, and spiky or bipo- 
lar), only diverging shocks can be evolved from such a sys- 
tem modulated by the Mach number M. 



parallel electric field decreases away from the center of 
the flow (see the amplitude decrease of Ei from Figl3] 
to FigU). In this way, we anticipate that an alternative 
solution would yield the classical inverted-^ structures, 
with scale length comparable to the ion gyroradius. 

Acknowledgements This work is funded by Visit- 
ing Fellowships in Canadian Government Laboratories 
Program, Natural Sciences and Engineering Research 
Council of Canada. The author thanks Senior Scientist, 
W. Liu, Canadian Space Agency, for valuable comments 
and suggestions. 
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Fig. 7 Nonlinear structures detected by the four Cluster 
satellites at four different locations and an ~0.03 s interval 
from each satellite, with amplitudes often reaching 500-750 
mV/m, the largest amplitudes measured in the outer mag- 
netosphere. Notice that SCI panel illustrates converging 
shocks, while others record diverging shocks. Adapted from 
Fig. 11 of Ref.[39]. 
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